Spin-1/2 frustrated antiferromagnet on a spatially anisotopic square lattice: 

contribution of exact diagonalizations. 
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The phase diagram of a spin- 1/2 J — J' — J2 model is inves- 
tigated by means of exact diagonalizations on finite samples. 
This model is a generalization of the Ji — J2 model on the 
square lattice with two different nearest-neighbor couplings 
J, J' and may be also viewed as an array of coupled Heisen- 
berg chains. The results suggest that the resonnating valence 
bond state predicted by Nersesyan and Tsvelik [Phys. Rev. B 
67, 024422 (2003)] for J2 = 0.5J' < J is realized and extends 
beyond the limit of small interchain coupling along a curve 
nearly coincident with the line where the energy per spin is 
maximum. This line is likely bordered on both sides by a 
columnar dimer long range order. This columnar order could 
extends for J' — > J which correspond to the Ji — J2 model. 
PACS numbers: 75.10.Jm: 75.50.Ee; 75.40.-s 



I. INTRODUCTION 

The interplay between frustration and quantum fluc- 
tuations in antiferromagnetic (AF) spin-1/2 systems in 
two dimensions (D) has attracted much experimental and 
theoretical attention in past years. In particular, many 
studies have been devoted to identify the possible phases 
which may appear at T = after the destabilization of a 
coUinear AF Neel phase by quantum fluctuations when 
increasing frustration in SU (2) invariant spin models. 
Large- approaches predict then a valence bond crys- 
tal (VBC) phase, with long range order (LRO) of sin- 
glet entities (dimers, plaquettes of resonnating dimers...) 
and thus breaking translational symmetry and/or others 
space symmetries, if the transition out of the Neel phase 
is continuous [1-3]. Numerical results obtained from ex- 
act diagonalization (ED) studies support this prediction 
for several models: the Ji — J2 — J3 model on the hon- 
eycomb lattice [4] and some models where the spins sit 
on a square lattice such as the crossed chain model [5,6] 
and the quadrumerized Shastry-Sutherland model [7] . 

Yet, for the much studied Ji — J2 model on the square 
lattice with a next-nearest-neighbor frustrating coupling 
J2 in addition to the nearest-neighbor coupling Ji (see 
Fig. 1 where the Ji — J2 model corresponds to J' = J = 
Ji), the validity of this scenario remain unclear [8-17]. Its 
classical ground-state has (tt, tt) Neel LRO for J2 < 0.5 Ji 
and is the degenerate manifold corresponding to two de- 
coupled Neel states on interpenetrating square lattices 
for J2 > 0.5Ji. There is a large consensus that (tTjTt) 
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FIG. 1. Exchange interaction pattern in the J — J' — J2 
model. The spins sit at the vertices shown by bullets, the 
coupling constant is J between nearest-neighbor pairs on 
hozirontal bonds (full lines), J' on vertical bonds (short 
dashed lines) and J2 on the diagonal bonds (long dashed 
lines) . The Ji — J2 model corresponds to J' = J = Ji . 



Neel LRO (see Fig. 2) survives to quantum fluctuations 
for J2 ^ O.4J1 and quantum fluctuations select either 
a (tt, 0) or a (0, tt) Neel state in the manifold of classical 
ground-states for J2 > 0.65Ji, whereas, in the intermedi- 
ate range 0.4 Ji < J2 ^ 0.65 Ji between these two gapless 
states, there is a phase without Neel LRO gapped to spin 
excitations. The transition out of the (tt, tt) Neel phase 
is second order while the transition out of the (tt, 0) Neel 
phase has been predicted to be of flrst order [10] but may 
be close to second order [14]. 

Whether this intermediate phase has VBC LRO or is 
some resonnating valence bond (RVB) phase which do 
not break any space symmetries, i.e. a spin-liquid, is 
debated. Earliest investigations, from ED calculations 
on samples up to A^ = 36 spins, have shown enhanced 
columnar dimer-dimer correlations in this intermediate 
region, but the finite-size scaling of the order param- 
eter for the columnar VBC LRO of dimers, shown in 
Fig. 3(a), were not conclusive [8,9]. Dimer series expan- 
sions, however, were in favor of VBC LRO: earliest stud- 
ies favored columnar LRO [10,11] while latter calcula- 
tions have lead to propose the occurence of two different 
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FIG. 2. Spins orders in the classical (tt, n), (tt, 0) and (0, tt) 
Neel states. 
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FIG. 3. Patterns of VBC orders proposed for Ji — J2 model: 

(a) columnar dimer, (b) plaquette, (c) columnar dimer with 
plaquette type modulation. Full lines connect spins bound 
into singlet: two spins in (a), four spins in (b) and (c). In 

(b) the strengh of spin-spin correlations on vertical and hori- 
zontal links of a plaquette i,j, k, I are the same but differs in 

(c) . There are four degenerate ground-states for orders (a), 
(b) and eight for the order (c) in the Ji — J2 model. These are 
obtained by one step translations in the horizontal and verti- 
cal directions and/or a 7r/2 rotation around a lattice site. 



VBC phases: the columnar dimer phase for J2/J1 < 0.5 
and a columnar dimer phase with plaquette type modula- 
tion, as shown in Fig. 3(c), for 0.5 > J2/J1, separated by 
a second order transition [14]. On the other hand, stud- 
ies based on the Green function Monte Carlo method, for 
systems of sizes larger than those which may be investi- 
gated by ED calculations, first favored either a plaquette 
LRO [12], as shown in Fig. 3(b), or a columnar dimer 
LRO with plaquette type modulation [13], similar to one 
of Ref [14], for J2/J1 = 0.5. Besides, it was concluded 
in Ref [12] that the finite size scaling of the susceptibil- 
ity to columnar dimer LRO from both Monte Carlo and 
ED calculations would exclude columnar dimer LRO at 
J2 — O.5J1. However, latter calculations [15], by the 
same group which favored the plaquette LRO [12], have 
lead to conclude in favor of the occurence of a spin-liquid 
phase for 0.4Ji < J2 < 0.5 Ji whereas it was pointed out 
that the ocurrence of the plaquette LRO will be unlikely 
in view of the spatial symmetries of the low lying singlets 
in the ED spectra. In their latest study these authors 
suggest either a spin-liquid or, possibly, the columnar 
state but with a weak order parameter [17] Yet, as the 
Monte Carlo calculations and the series expansions start 



from some trial state which may bias the results and in- 
volve approximations which are not fully under control, 
while ED calculations which were performed on samples 
having different lattice symmetries and display irregular 
finite size scaling, the nature of the ground-state in the 
intermediate region remain an open question. 

In this paper, motivated by a recent paper of Nersesyan 
and Tsvelik [18], refered as (NT) in the following, we at- 
tempt to shed light on this issue from ED calculations on 
a spatially anisotropic version of the Ji — J2 model with 
two different nearest neighbor couplings. The Hamilto- 
nian of this J — J' — J2 model (named the confederate 
flag model by NT) reads: 

+ 

)] (1-1) 

which may also be viewed as an array of M spin chains of 
lenght L, if J' < J. The exchange J (see Fig. 1) couples 
first-neighbor (horizontal) pairs of spins along the chains, 
exchange J' first-neighbor spins in the (vertical) trans- 
verse direction and exchange J2 second-neighbor pairs on 
the diagonals of the square plaquettes. All exchanges are 
assumed positive, describing AF couplings. For J' = J 
one recovers the Ji — J2 model. Mostly interested with 
the limit M, L — > 00, we shall also focus on the case 
J' < J. 

The classical ground-state of the J — ,/' — J2 model is 
rather similar to the one of the Ji — J2 model. When 
J' < J, one has {n,n) Neel LRO if J2 < 0.5 J', (tt, 0) 
Neel LRO if J2 > 0.5 J' and the (horizontal) J chains are 
decoupled if J2 — 0.5 J'. For J' > J, one has (tTjTt) Neel 
LRO if J2 < 0.5 J, (0,7r) Neel LRO if J2 > 0.5 J whereas 
the (vertical) J' chains are decoupled if J2 — 0.5 J. 

This Hamiltonian has been much studied in past years 
for a finite number of chains, the so-called M-leg ladders, 
generally in the case of open boundary conditions in the 
transverse direction (for M > 2 since if M = 2 open and 
periodic boundary conditions are equivalent). 

Numerical studies [27-29] based on the density matrix 
renormalization group (DMRG) method or Monte Carlo 
approaches have shown that the M-leg unfrustrated lad- 
ders (J2 = 0) display a behavior analogeous to a M/2- 
spin chain [26], being fully gapped if M is even and gap- 
less if M is odd. So (tt, tt) Neel LRO only occurs in the 
limit M — > cxo. 

The effect of frustration has been especially investi- 
gated theoretically and numerically for the 2-leg ladder 
(see Ref [21-24] and references therein). These studies 
have concluded that one has two phases separated by a 
transition line which has been found to coincide with the 
line of maximum frustration where the energy per spin 
reaches its maximum, noted J2^{J') in the following for 
all values of the number of chains. This line starts as 
J^'(J') - 0.5 J' for J' and bend shghtly with in- 
creasing J' so that J^(J') - 0.6 at J' = 1 [21,22] (as 
the dashed hue in Fig. 4). If J2 > J2^{J') one has the 
'Haldane phase', where the two spins on a rung tend to 
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be in a triplet state, so named as it contains the point 
J2 = J, J' = whose low-energy speetruni is similar to 
that of a spin-1 Hcisenbcrg chain. If J2 < J™ {J') on has 
the so called 'singlet phase', so called as it contains the 
case J' ^ J with J2 = 0, where the ground-state con- 
sists of singlets on each rung. As shown by White [31] 
the 'singlet phase' is also an Haldane phase with diago- 
nally situated next-nearest neigbor spins coupling to form 
an effective S = I. Both phases have a topological or- 
der, breaking a hidden Z2 x Z2 symmetry, which can be 
measured by non-local 'string order' parameters, similar 
to the string order parameter of the S = 1 Heisenberg 
chain. The two phases are fully gapped like the S = 1 
chain [26], and have a non degenerate singlet ground- 
state. These ground-states become degenerate on the 
transition line J™ (J'). The nature of the transition has 
been debated. In a DMRG study [22] it was predicted 
that the transition line would be gapless for J2 lower than 
~ 0.287 J. But more recent DMRG calculations [23] have 
concluded that a gap subsits between the two-fold degen- 
erate ground-state and the excited states for all J2 > 0. 
So the transition is 1th order on the whole line J2'{J') in 
agreement with bosonization results in the limit of van- 
ishing interchain coupling [24]. The elementary excita- 
tions then consist of gapped deconfined spin-1/2 spinous 
which may be viewed as kinks interpolating between the 
two ground-states [24] . Note that these results imply that 
the classical behavior of independent chains on the line 
J2 = 0.5 J' is then destabilized by quantum fluctuations. 
The single Heisenberg chain is gapless. The spectrum of 
two independent chains too. The existence of a finite gap 
excludes such a behavior. 

A DMRG study of the three-leg frustrated ladder with 
open boundary conditions in the transverse direction [25] 
has shown that this system exhibit two phases separated 
by a transition line also coincident with the line of maxi- 
mum frustration J^(J') and which is a curve quasi iden- 
tical to the one found for the 2-leg ladders. The phases 
are analogeous to those of the 2-leg ladders. The small 
J2 phase, which is the phase of the unfrustrated ladder, 
show a tendency of the three spins on a vertical line to 
pair in a state of minimum spin whereas they tend to pair 
in a state of maximum spin in the large J2 phase, which 
is equivalent to the spin-3/2 chain. On the transition line 
a bosonization study at weak interchain coupling [32] has 
predicted that the ground-state is a chiral spin-liquid. 

Recently NT, using an approach based on the 
bosonization method, investigated the crossover from fi- 
nite M to M ^ 00 in the J — J' — J2 model on the 
line J2 = G.5J', in the limit of small interchain cou- 
pling (J', J2 <C J) and predicted that the classical be- 
havior of independent chains is unstable to a special 
RVB (spin-liquid) state which may be a realization of 
the chiral 7r-flux RVB state [18-20]. Its ground-state 
has a 2*^~^ degeneracy, for M even and transverse pe- 
riodic boundary conditions, associated to the breaking 
of a local Z2 symmetry, present in the bosonised ver- 
sion of the model, corresponding to the invariance under 



independent translations by one lattice spacing along in- 
dividual chains. Each ground-state may be caracterized 
by different values of a set of non-local order parame- 
ters. These non-local order parameters correlate spins 
on two adjacent chains whereas different pairs of chains 
are uncorrelated. The elementary excitations consist of 
deconfined spin-1/2 spinous interpolating between differ- 
ent ground-states. These excitations are either gapless, 
as predicted for some values of the number of chains, 
such as M = 4,6, 12, or gapped, as for the 2-leg ladder 
(M = 2). 
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FIG. 4. Proposed phase diagram of the J — J' — Ji model 
for L — > 00 and M— >oo(J=l). The solid lines indicate 

the approximate location of the boundaries of the (tt, tt) and 
(tt, 0) Neel phases, noted J^{J') and J^{^J') respectivelly. 
The dashed line in the intermediate region between the Neel 
phases is the approximate location of the maxima J™ {J') of 
the energy per spin eo along which the present ED results 
suggest that the RVB state predicted by NT may be realized 
up to large values of J'. In the intermediate region, out- 
side the curve J^, the ground-state is predicted to display 
columnar dimer order. The Neel phases only appear for both 
L, M — > 00. For finite M and L —f 00 one has gapped phases. 
The position of the curve J™ (J') is quasi-independent of the 
number M of chains down to M = 2, which correspond to the 
2-leg ladder. 

Two questions then arise. First, does the state pre- 
dicted by NT extends outside the range of small inter- 
chain coupling. Second, what is the nature of the behav- 
ior in its vicinity, in particular for the case of an infinite 
number of chains. Below we report results of ED calcu- 
lations oi N = M X L spins for < 36 to shed light on 
these questions and more generally the phase diagram of 
the J — J' — J2 model. We shall consider only the case 
of an even number of chains. The ED calculations have 
been carried out on samples oi N = M x L spins in M 
chains of lenghts L with periodic boundary conditions 
both along and perpendicular to the chains: samples of 
4 chains with N = 16, 24, 32, 36 spins and 6 chains with 
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N = 24, 36 spins. Additional calculations were also per- 
formed for the M = 2 ladder for < 32 in order to 
compare with the M = 4, 6 results. Contrary to the 
samples considered in previous ED calculations for the 
Ji — J2 model, all samples have translation vectors par- 
allel to the basis vectors of the square lattice which leads 
to a regular evolution of the properties with increasing L 
and M. 

Examination of the ED spectra indicates that the 
model will display two Neel phases, with magnetic wave- 
vectors (tt, tt) and (tt, 0), separated by a magnetically dis- 
ordered intermediate region, in the limit A — > cx) with 
both L — > cx) and M 00, as shown in Fig. 4. At 
variance with a study based on the DRMG method [33] , 
which appeared very recently on a preprint server, our 
results provide indications that the classical behavior of 
independent chain is destabilized by quantum fluctations 
and the RVB state of NT may occur along a curve co- 
incident with the line of maximum frustration J2^{J'), 
located in this intermediate region, up to large interchain 
coupling. They also indicate that the intermediate region 
of the M 00 phase diagram (see Fig. 4), on both sides 
of this line, has columnar VBC LRO, which may extend 
up to J' — > J and could already appear for M even > 4. 

In Sec. II we describe the evolution of the energy per 
spin vs J2 and J' and the location of the curve J^{J'). 
The properties of the model in the regions of (tt, tt) Neel 
LRO, (tt, 0) Neel LRO and in the intermediate region are 
presented in Sec. Ill, IV and V, respectively. Sec. VI 
gives a summary of the our results and discuss briefly the 
possible behavior of the model for ferromagnetic J2 and 
J'. 
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FIG. 5. Ground-state energies per spin eo = Eo/N vs J2 
for the N = 24: sample with M = 4: [4 chains). The lines are 
guides to the eyes. 
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FIG. 6. Same as Fig. 5 for the iV = 24 sample with M = 2 
(2-leg ladder). 



II. ENERGY RESULTS 

The ground-state energies per spin bq = Eq/N are 
ploted, as a function of J2, for different values of J' < 1 
( J = 1 in the following) for the AI — 4 chains, A^ = 24, 
sample in Fig. 5. For comparison, the results for the 
A^ = 24 2-leg ladders are indicated in Fig. 6. Fig. 7, 
displays eo vs J2 at J' = 0.8 for the different ill = 4 
and M = 6 samples. As shown in Fig. 5 and Fig. 7, 
the point of "maximum frustration" where eo reaches 
its maximum, J™ (J') occurs for a value of J2 slightly 
larger than 0.5J', as found for the two and three-leg lad- 
ders [21,22,25]. 
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FIG. 7. Ground-state energies per spin eo = Eo/N vs J2 
at J' =0.8 for different samples of M chains of lenghts L. 
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As shown in Fig. 7, J™ slightly increases with L, 
but is nearly independent of the number Af of chains. 
This indicates that the L ^ oo curves J™ (J') will be 
quasi-independent of M, as previously conjectured from 
a comparison of the results for the two and three-leg lad- 
ders [25] . The maxima of eo (for a given sample) decrease 
with increasing J' but remain rather close to the value 
for decoupled chains at J' = at least if J' < 0.9. 

In Fig. 5 and Fig. 6 one sees, in the range 0.9 < J' < 1 
a drop of eo at J™ from a value close to the one for decou- 
pled J chains of lenghts i to a value closer to the one for 
independent J' chains of lenghts M. This corresponds to 
the crossover to the situation at J' > 1 where the finite 
samples M x L are best viewed as L coupled chains of 
finite lenghts M. 

III. (tt, tt) NEEL long range ORDER 

Increasing J2 from zero at given J' (< 1), one finds 
a range of parameters, which extends to a value we note 
J|i(J'), sHghtly smaller than 0.5J' (and thus < JTiJ')), 
where the spectra exhibit the features at finite size spe- 
cific of a system displaying coUinear (tt, tt) Neel LRO in 
the limit N —>■ (x with both L —>■ 00 and M —>■ 00. The 
Neel order breaks SU{2) and lattice spacial symmetries. 
Evidence of Neel LRO in the spectra is the presence of 
lowest eigenlevels in each spin S sector which form a set of 
states with energies E{S) - Eq S{S -I- l)/iV for i -> 00 
and M 00, the so-called quasi degenerated joint states 
(QDJS) [34,5] as shown in Fig. 8. These QDJS are 
well separated, at finite size, from the others lowest ex- 
cited states (magnons) and collapse on the ground-state 
faster with N than the magnons, enabling the breaking of 
SU{2) and lattice spacial symmetries. For coUinear LRO 
there is one QDJS for each S value. The QDJS specific 
of (tt, tt) Neel LRO consist of a state belonging to an ir- 
reducible representation (IR) with wave-vector k = if 
S is even or k = (tt, tt) if S is odd, both even in a spatial 
R{tt) rotation around a site and in the reflexion <j with 
respect to a J chain. 

The approach of ED calculations does not allow an 
accurate location of boundary of the Neel phase since 
it is limited to small samples. Nevertheless the exten- 
sion of the Neel phase can be approximatively estimated 
from ED calculations by monitoring the range of val- 
ues of J2 for which the QDJS remain well defined in 
the spectrum [5]. This gives J|"'^(J' = 1) « 0.4 for the 
Ji - J2 model, J|i(J' = 0.8) « 0.35. J|i(J') moves 
steadily closer to the line J2 = 0.5 J' as J' decreases. 
Besides, since quantum Monte Carlo calculations on the 
unfrustrated line of the phase diagram (J2 = 0) have 
indicated that (tt, tt) LRO appear as soon as J' > 
for M — *■ 00 [30], Neel LRO extend down to the point 
J' = 0, J2 = 0. Thus the curve Ji^iJ') is tangent to the 
line J2 = 0.5 J' for J' 0, moving away from this line 
as J' increases from zero, and is located below the line. 
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S(S+1) 

FIG. 8. Spectrum at J' = 1 and J2 = vs S{S + 1) for 
A'' = 36. The QDJS characteristic of coUinear Neel LRO at 
the bottom of each S sector are well aligned (dashed line) and 
clearly separated from the other lowest excitations (see text). 

as shown in Fig. 4. 

It is to be emphasized that Neel LRO appear only for 
both M — > 00 and L — » 00. There is no Neel LRO 
for a finite number M of chains. Indeed, as mentioned 
above, systems which consist of a finite number M oiS = 
1/2 Heisenberg chains are known to acquire a gap upon 
shifting a tranverse coupling J' if M is even, behaving 
analogeously to a spin S = M/2 chain, so these systems 
are in a Haldane phase [27-30]. 

The typical evolution of the spin-gap with M and 
L is illustrated in Fig. 9 where values of the A-'^ are ploted 
vs 1/7V for N = 16, 24, 32 samples of M = 4 chains and 
the TV = 24, 36 samples of M = 6 chains at J' = 0.8 and 
J2 = 0. The evolution for M = 4 of A^ vsl/N ^ 1/L 
shows a small but noticable upward curvature as found 
previously for the 2-leg and 4- leg ladders (see Fig.l of 
Ref. [27]). As found in studies of the 2,4,6 legs ladders, 
one may assume that the gap decreases monotonically 
as 1/L decreases with a curvature that remains positive 
and vary only smoothly. To estimate a lower bound on 
A^{L — > 00) we use a linear extrapolation vs 1/L, be- 
yond > 32 of a spline fit to the M = 4 three values of 
A^ (dashed line in Fig. 9). This leads to a lower bound 
of A^ for L 00 ^ 0.09. This lower bound indicates 
that A^{L — > 00) remains finite. Due to the neglect of 
the curvature beyond N > 32 the gap is slightly under- 
estimated. The true value of A^ for L ^ 00 for M = 4 
is likely somewhat closer to the value (~ 0.14) estimated 
for the 4-leg ladder with open boundary conditions in 
the transverse direction [28]. For M = 6, the present 
data, limited to the two values for L = 4, 6, show that 
A^ {L — > 00) decreases with increasing M, and remain fi- 
nite if M = 6. Its value for the 6-leg ladder is ~ 0.04 [28] 
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FIG. 9. Spin-gaps at J' = 0.8 and J2 = vs 1/iV for for 
— 16, 24, 32 samples of M = 4 chains (squares) and the 
N = 24, 36 samples of M = 6 chains (diamonds) . The dashed 
line is a spline fit to the M = 4 values for N <32 followed by 
linear extrapolation for N > 32. 



if J' = 0.8. The present data show that the gap is very 
small but do not allow an accurate estimation of a gap 
of this order of magnitude. 

The evolution of the spin-gap with N ^ 00 for M = L 
is different than for fixed M and L — > 00. For M = L, 
where ^ for A'' — > 00, the N = 16,36 values in 
Fig. 9 indicates an evolution of the spin-gap in the case 
M = L that decreases with a negative curvature vs 1/N 
as usually found in systems with Neel LRO [35]. 



IV. (tt, 0) NEEL LONG RANGE ORDER 

For J2 larger than a value J|^, larger but very close to 
J2", one similarly recognizes in the spectra, the features 
specific of collincar Neel LRO with a wave vector (tt, 0) 
or (0, tt) depending on value of J' /J. 

li J' = J = 1 and M ~ L these two Nccl orders are 
degenerate. The QDJS then include two states for each 
value of S: two with k = if S is even and two with 
k = (tt, 0) and k = (0, tt) if S is odd, all even in -R(7r) 
and a. 

The degeneracy is lifted \i J' ^ J. For J' < J and in 

the thermodynamic limit M = L ^ 00 the (tt, 0) LRO 
is favored whereas this would be (0, tt) LRO for J' > J. 
The QDJS of (tt, O) LRO are observed as the lowest states 
of the spectra as soon as J' < J for M = L. They include 
only one state for each value of S: a state with k = if 
S is evcin and a state with k = (tt, 0) if S is odd, both 
even in R{Tr) and a. 

Note, however, that at finite size, for an aspect ratio 
M/L < 1 of the samples, quantum fluctuations, which 



favor an antiferromagnetic allignement of the spins in the 
shortest direction, lead to spectra caracteristic of (0, tt) 
LRO even for J' < J as J' w J. The transition between 
the two kind of spectra then occurs in the range of values 
of J' where eo at J™ (J') drops to a value similar to the 
one of L decoupled chains (see Fig. 5). For M = 4 and 
L > M, it is only for J' < 0.9 that the QDJS of (tt, O) 
LRO becomes lower than the QDJS of (0,7r) LRO. 

The boundary of the (tt, 0) Neel phase can be also es- 
timated in the same way as for the (7r,7r) phase. For 
the Jl — J2 model it was found that Jo^ occurs at 
a value slightly larger than J2 
0.60 < JfiJ' = 1) < 0.70 



0.6, in the range 



tra for J' < 1 indicates that move closer to 



Examination of the spec- 
J' 

decreases with J^^ « 0.5J' for J' 0. As at J2 = 0, for 
Heisenberg chain coupled by J' interactions, it is likely 
that, at J' = 0, the ID behavior of the single Heisen- 
berg chain is unstable to J2 interchain coupling leading 
to Neel LRO as soon as J2 > 0. 

Here too, the Neel LRO appears only in the limit of 
an infinite number of chain. The spin-gap remains finite 
at M fixed for L — > 00. The ground-state is a non de- 
generate singlet. One has Haldane-likc phases. Instead 
of the (tt, tt) Neel LRO one has the same phase as the 
phase of the M-lcg (unfrustrated) ladder. In place of the 
(tt, 0) Neel LRO one has a phase, that is probably closely 
similar to the one of the spin-M/2 chain as found for the 
2-leg and 3-leg ladders. These phases are analogues of 
the 'singlet' and 'Haldane' phases of 2-leg ladder. Like 
the 2-leg ladder they have likely an hidden topological 
LRO. This and their eventual relation to the chain of in- 
teger spins S* > 1 for which string order parameters have 
been recently studied [36] deserves further investigation. 



V. INTERMEDIATE REGION 

Between these two regions of Neel behavior, there is an 
intermediate region for the range of parameters ^ {J') < 
J2 < J2^{J') where the spectra indicate an absence of 
Neel LRO when L 00 and M — > cxd. Leaving the region 
of (tt, tt) Neel LRO by increasing J2 at constant J', the 
evolution of the energies of the lowest eigenstates in each 
S sector changes to E{S) — Eq ^ S for small values of S 
(see Fig. 10). This feature is an indication that the spin- 
gap opens for L ^ cx) and M — > 00 [37] and one enters 
a magnetically disordered region. At the same time one 
observes a lowering of some singlet states. This linear 
behavior reaches its maximum extension for a value of 
J2 that appears to coincide with J™ where it extents up 
to S* = M/2 if J' < 0.9 as shown in Fig. 11 for = 4 at 
J' = 0.8. Then, beyond J™, the evolution of the lowest 
eigenenergies returns fast to -^(5*) — Eq ^ S{S -I- 1) in the 
(tt, 0) Neel region which starts at J^^, close to J2'. 

The quantum numbers of the first excitations in each 
spin sector and the signs of the spin-spin correlations 
remain the same as in the (tt, tt) Neel region for J2 < J™ 
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FIG. 10. Spectrum of eigenenergies vs total spin S for the 
M = 6 sample of A'" = 36 spins at J' = 0.6 and J2 = 0.3. 
These values of J' and J2 correspond to a point in the in- 
termediate region (shghtly) below the line J2"(J') in Fig. 4 
{J^fyJ' = 0.6) PS 0.32). The dashed hne, fitted to the lowest 
5 = 1 and S = states, shows that the energy of the lowest 
states in each spin sector increase nearly as ~ S at small S 
values. 
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FIG. 11. Spectrum at J' = 0.8 and J2 = 0.45 for the 
N = 32 sample of M = 4 chains 
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FIG. 12. Spectrum of decoupled chains for the A'^ = 32 
sample of M = 4 chains (J' = and J2 = 0). Note that the 
excited levels have a large degeneracy. 



and are the same as in (tt. 0) region for J2 > J^. The 
curve J^^i J') is the analogue of the transition hne of the 
classical model J2 — 0.5J'. On this line the classical 
ground-state consist of decoupled chains. Then for J2 w 
J^{J') one may notice features, in the spectra and the 
spin-spin correlations (observed for M = 2,4,6), that 
present certain similarities with those of M decoupled 
chains if J' < 0.9 (and those of L independent chains 
when L> M and J' > 0.9). 

As indicated above, one has E{S) — E'o ~ "S* up to 
S = M/2, as for decoupled chains (see Fig. 12), which is 
what would result for independent magnetic excitations 
on the chains. Moreover, the very lowest eigenstates in 
each S sector belong to the same IR as those of a sys- 
tem of decoupled chains. For instance, the first triplets 
excitations is a set of M states with all wave-vectors of 
type (tt. ky) on the side the Brillouin zone which become 
quasi degenerate at J2 « J™ as shown in Fig. 13 for the 
M = 4, TV = 24 sample at J' = 0.8. Similarly, first ex- 
citations at S' > 1 and S < M/2 have same quantum 
numbers (wave vector, characters in point-group symme- 
tries) as those of decoupled chains which result from the 
combination of first triplets excitations. The energies per 
spin (see Fig. 7) and the spin gaps at finite size are also 
quasi independent of the number of chains. 

In addition the spin excitations appear quasi gap- 
less. Although the finite-size spin-gaps are maximum for 
J2 ~ J™, as shown in Fig. 13, the analysis of the evo- 
lution of the spin-gap of the M = 4 samples with 
1/L indicates that its value at L — > 00 is minimum for 
J2 w as for the 2-leg ladder [22,23]. If J' < 0.6 the 
scaling behavior of is quite similar to the one for the 
single chain (its evolution with 1/L is quasi linear with a 
very small negative curvature) and the extrapolated value 
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J2 J2 



FIG. 13. Gaps to the 1th triplet (AS = 1) vs J2 for AT = 24 
samples with M = 4 (4 chains) and wave vectors k = {n,ky) 
on the side of the Brillouin zone if J' — 0, 0.4, 0.8. The lines 
are guides for the eye. 

of is « 0. At J' > 0.8 nevertheless the evolution of 
with 1/L turns to have a positive curvature and A^ 
extrapolates to a finite (but small) value. This may be 
consistent with a spin-gap opening exponentially with 
increasing J' along J™ (J') as predicted by NT at weak 
interchain coupling along J2 = 0.5J'. However the spin- 
gap on J™ (J') appears to remain very small up to large 
interchain coupling for M < 4 as found for AI — 2 [23] . 

An apparent decoupling of the chains is also manifest 
in the spin-spin correlation between two spins at distance 
(/,m): 

s(;,m) =< So, o-S;,m > . (5.1) 

Values of s{l,m) for the = 32 (M=4) sample are dis- 
played vs J2 at J' = 0.8 in Fig. 14 for pair of spins on the 
same (m = 0) or different (m = 1, 2) chains at short dis- 
tances along the chains {I < 2). Since s{l,m) decrease in 
magnitude at given m with increasing distance I along the 
chains these are the largest values of s{l,m). J™ « 0.45 
when J' = 0.8 and N = 32. As shown in Fig. 14 the 
inter-chain spin-spin correlations become very small at 
J2 = 0.45: if m ^ 0, s{l,m) w except between first 
neighbor spins which changes of sign at a slightly larger 
J2 value. 

Others features in the spectra nevertheless testify of 
interactions between the chains when J2 « JT{J') f^^' 
M < 4 as found for M = 2. Starting from situation of de- 
coupled chains at J' = and increasing J' along J™ (J') 
one sees a modification of the structure of the eigenlevels 
above the most lowest eigenstates in each S sector. Levels 
which were degenerate for uncoupled chains move away 
from one-another. In particular, there is a number of sin- 
glet states which separate from the others and drift to- 



FIG. 14. Spin-spin correlations s{l,m) (Eq 5.1) between 
spins at distance (/, m) vs J2 for = 32, M = 4, J' = 0.80. 
The lines are guides for the eye. 

ward the ground state (from an energy w twice the spin- 
gap at J' = to an energy ^ to the spin-gap at J' — 0.8 
as shown in Fig. 11 for the iV = 32(M = 4 x L = 8) 
sample). This could indicate a behavior at J2 ~ J™ 
that differs from the one of independent chains and raises 
the question wether it corresponds to the RVB state pre- 
dicted by NT in the limit of small interchain coupling for 
J2 = 0.5J'. 

We investigate below successivelly the three cases J2 < 
J™, J2 ss and J™ > J2 in more detail. In order to 
be with an intermediate region sufficiently wide and to 
observe effects of interactions between the chains on small 
samples we shall mainly display results for a rather large 
J' = 0.8 value and discuss their evolution for smaller and 
larger values of J'. 

A. J2 < 

For J2 < J™ and not to close to J™, the lowest excited 
singlet states, noted below |1 >, |2 >, |3 > are: |1 > in 
the trivial IR (quantum numbers: k = 0, i?(7r) = 1,(t = 
1), |2 > [k = (7r,0),i?(7r) = -l,cr = 1] and |3 > [k = 
(0,7r),i?(7r) = -l,cr = -1]. 

Fig. 15 shows the evolution of the singlet-gaps A^ from 
the ground-state |0 > to the states \i > vs 1/L at J' — 0.8 
and J2 = 0.4 for the M = 4, = 16, 24, 32 samples. The 
gaps appear to decrease monotonically. They evolve with 
an upward curvature vs 1/L. As above we assume that 
the evolution of the gaps will remain similar at larger 
sizes as found for the 2-leg ladder. Lower bounds for 
the gaps are then obtained using linear extrapolations 
for A^ > 32. The lower bounds are positive values for A" 
and A3. This point to finite gaps A]', A3. On the other 
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hand, the extrapolation suggests a possible vanishing of 
A2 for L ^ 00. A similar analysis of these gaps for 
the M ^ 6, N = 24,36 samples also point to the same 
conclusions for L — > 00 if M = 6. 

The vanishing of A2{L — > 00) for M < 4 may be 
questioned since the 2-leg ladder (M = 2) is known to 
be fully gapped [23]. The same analysis indeed shows 
that A2(L —>■ 00) is finite and rather large if M = 2 at 
J' = 0.8 and J2 = 0.4. From this, it could be argued that 
A2(L 00) decreases with increasing values of M but 
remains finite for M finite. Our data, limited to small 
sizes cannot exclude this possibility. However, A2 most 
likely vanish on J™ for M > 4 (sec Sec. V B) and then, 
the range of values of J2 where A2 appear to vanish, 
starts at J2 « 0.35 and extends beyond J™ till J2 « 0.5 
(see below) if M = 4. This interval is large which sup- 
port a vanishing of A2 on an extended range of values of 
J2. An other possibility, is that the value of J2 at given 
J' beyond which A2(L — * 00) would vanish vary with M 
and drift toward J™ as M decreases. ED calculations 
indeed point out that A2(L — > 00) drops to a minimum 
at J2 « J™ for M = 2 where this gap is quasi- vanishing. 

On the other hand, an analysis of the evolution of the 
spin-gap A^ with 1/L show that it remain finite if i — > 00 
for M = 4 and J2 < J^". Besides the fact that E{S) 
neither evolve as E{S)-Eo ~ S{S+1) or as E{S)-Eo ~ 
S indicate that the spin-gap is finite for L — » 00 and 
M ^ 00. 

All together, the ED results suggest the occurence of 
a (77,0) VBC LRO at J' = 0.8 for J2 > 0.35, breaking 
translational symmetry in the horizontal direction with a 
twice degenerate ground-state, gapped to others singlet 
and magnetic excitations, possibly as soon as M = 4 
when M is even, and most likely when M 00. 

The patterns of dimcr-dimer correlations 

D{i,j;k,l) =< (Si.S,)(Sfc.Si) > - < Si.Sj >< Sk.Si > 

(5.2) 

between valence bonds on pairs sites and {k,l) are 
displayed in Fig. 16 for the = 32 sample at J' = 0.8 
and J2 = 0.4. Fig. 16(a) is consistent with a colum- 
nar (tt, 0) VBC LRO although the correlations between 
horizontal dimers on different chains are rather small. 
Fig. 16(b) shows nonetheless correlations between verti- 
cal dimers that are not smaller than those between hori- 
zontal dimers and that could be consistent with a colum- 
nar (0, tt) VBC LRO. This may be related to the fact that 
A2 at L = 8 is still higher than the singlet gaps A°,A3. It 
is the extrapolation of A", A3 at L —>■ 00 to finite values 
that exclude a (0, tt) VBC LRO for J' = 0.8. In Fig. 16 
we see that the dimer correlations are strongest between 
pair of bonds on adjacent chains and much weaker other- 
wise. These features reinforce with increasing J2 and will 
be maximum at J2 ~ J™. Simultaneously, the extrapo- 
lated values of AjjAg at L — ^ 00 decrease toward zero. If 
J' < 0.9, the same features for the singlet gaps and the 
dimer correlations remain in the intermediate region for 
J2 < J,". 
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FIG. 15. Singlet-gaps A° at J' = 0.8 and J2 = 0.4 vs 1/L 
for for iV = 16, 24, 32 samples of M = 4 chains. The dashed 
lines are spline fits to the data for A'^ < 32 followed by linear 
extrapolation for JV > 32. 
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FIG. 16. Dimcr-dimcr correlation function D{i,j;k,l) 
(Eq 5.2) of a dimer on a reference pair of sites with a 

dimer on a pair of sites {k, I) for A'^ = 32, J' = 0.80, J2 = 0.4. 
The reference bond: (1,5) in (a), (1,4) in (6), is represented 
by a double line. A solid (dashed) line means a positive (neg- 
ative) value of D. The thickness of the line is proportional to 
the magnitude of D. 
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Interchanging the role of J and J' in the preceeding 
analysis points out that a (0, tt) VBC could occur for 
finite i if M ^ oo at sufficiently large J' . On the other 
hand, VBC LRO for finite M, in the limit L oo, may 
be excluded at large J'. As mentioned in Sec. II, there is 
a crossover in the range 0.9 < J' < 1 to a situation where 
the systems, for an aspect ratio M/L < 1, may be best 
viewed as an array of L chains in the vertical directions. 
For larger J' the evolutions of the above singlet gaps with 
1/1/ show that A2 opens and all gaps are finite at L — > cxd 
for M = 4, 6. The (tt, 0) VBC LRO then disappear. In 
particular, there is no VBC LRO for M finite if J' = 
J = I. The gaps nevertheless decrease at fixed L for 
increasing M and the value of J' where the crossover 
occurs for M/ L <1 shifts with increasing values of M to 
J' = J for M/L 1. Thus (tt, 0) VBC LRO may extend 
to J' ^ J if L, M ^ 00 where it would be degenerate 
with (0, tt) VBC LRO. 

However an analysis of the evolution of the gaps A°, 
A", A3 with the size of the samples one may study by 
the ED method do not allow to conclude that these gaps 
vanish at J' = J = 1 for A'' ^ 00, as required for the 
four-fold degenerate columnar order corresponding to de- 
generate (tt, 0) and (0,7r) orders. 

It could be possible that one has (tt, 0) LRO for J' < 1 
and (0, tt) LRO for J' > 1 but that VBC LRO vanishes 
on the line J' = J. However, the transition out of (tt, 0) 
VBC phase might be first order: the gap A2 appears 
to jump to a large value going through the transition 
by increasing J'. The order parameter would then be 
finite on the transition line and, if this transition line 
extends to the line J' = J for L = M ^ 00, it would 
imply that the four-fold columnar order occurs on the 
line J' = J where the J — J' — J2 model reduces to the 
Ji — J2 model. A picture of the dimer correlations similar 
to the one observed above for at J' ~ 0.8 and J2 = 0.4 
with strongly correlated pair of chains has been found 
in dimer series expansions approaches for the Ji — J2 
model when 0.4 < J2 ^ 0.5 [10]. The pattern of dimer 
correlations in one of the four-fold degenerate columnar 
states would then be reminiscent of the one off the line 
J' = J. 

The locations of the end of the (tt, tt) Neel phase (for 
M 00) and of the beginning of the (tt, 0) VBC phase 
may differs [16]. They arc difficult to estimate accuratclly 
from ED calculations with the present sizes but are prob- 
ably very close or coincident as conjectured for the Ji — J2 
model [16]. 

B. J2 fs J™ 

Jjp is « 0.45 at J' = 0.8. At J' = 0.8 and J2 = 0.45, 
the evolution of the gaps to the first triplet for M = 4 
with 1/L (not shown) indicates that the spin-gap remains 
finite L 00 although it has become much smaller than 
at J2 < J|". The new feature is the appearence of other 



singlet states in the low part of the spectrum which, at 
first sight, appear to form a set of states separated from 
the singlet contimuim (see Fig. 11). This set includes the 
lowest states with wave vector k = (0, ky) [ky = 27rm/M 
m e [M/2 — l,...,M/2]) even under a: the two states 
noted |4 > [k = (0,±7r/2), cr = 1] if M = 4 or the 
states [k = (0,±7r/3), a = 1], [k = (0,±27r/3), ct = 1] 
if M = 6, then noted |4 > and |4' >. If M = 4, there 
also the 2th excited state in the trivial IR, noted |5 >, 
just above the states |1 >, |2 >, |3 >, |4 > whereas if 
M = 6, there are also the lowest states [k = (tt, ±7r/3), 
a = —1], then noted |6 >. Next if M = 4, one has just 
below the singlet continuum a state, noted |6 >, with 
[k = (tt, ±7r/2), cr = —1], whereas if M = 6, one finds, 
adjacent to the singlet continuum the 2th excited state in 
the trivial IR (noted |5 >) followed by states with wave 
vectors k = (0, ky). 

This low energy spectrum is different from the one for 
independent chains to which it may be compared. In 
the case of decoupled chains the lowest singlets above 
the ground-state (see Fig. 12) is the degenerate set Si of 
M states with wave vectors k = (tt, ky) {ky = 2-Km/M 
m e [M/2 - 1, M/2]) which consists of |2 >, |6 > and 
the state [k = (tt , tt) , R{tt) — — 1,(t 1], noted |8 > if 
AI = 4 whereas it includes |2 >, |6 >, |8 > and the state 
[k = (7r,±27r/3), a = -1] noted |8' > if M = 6. These 
states correspond to 5 = lowest excitations along indi- 
vidual chains (which wave vector differs from the one of 
the ground-state by Afc = tt). Just above this set, one 
has the set 52 of degenerate states with wave vector 
k = (0, ky) which consists of |1 >, |3 >, |4 >, |5 > and the 
state [k = (0, tt), i?(7r) = 1, a = 1], noted |7 > if M = 4 
(6 states), whereas if M = 6, it includes |1 >, |3 > |4 >, 
|4' >, 1 5 >, 1 8 >, a 3th state in the trivial IR, a second 
state in the same IR as |4 > and two states in the same 
IR as |4' > (15 states). These states correspond to the 
combination of two S = 1 lowest excitations on different 
chains in a total singlet state (the lowest S = 1 excita- 
tions of the single chain has a wave vector which differs 
from the one of the ground-state by Ak = tt). Third, on 
has a set ^3 of M states with wave vectors k = {-K,ky). 
These states are: [k = (tt, 0), i?(7r) = — l,cr = — 1], 
[k = (7r,7r/2),a = 1], [k = (tt, tt), i?(7r) = -l,a = -1] 
if M = 4. They correspond to the combination of three 
lowest 5=1 excitations on different chains. Above 
appear sets of states with wave vectors inside the Bril- 
loun zone which involve higher excitations on individual 
chains at A/c ^ tt and appear to form a continuum. At 
small values of L, the sets of states <Si, 52 remain well 
separated from this continuum which starts at 53. Yet, 
in the limit L — *• 00 the gaps to the states of S\, S2, 
S3, vanishes as the gaps to the lowest 5 = 0, 1 excited 
states of a single chain. The excited states form a gapless 
continuum adjacent to the ground-state. 

On the line of maximum frustration, the degeneracy 
in the sets Si and 1S2 appears to be lifted by interchain 
couplings. Fig. 17 shows the variation of the gaps A° for 
the N = 32 sample vs J' for values of J2 « Jr(^')- The 
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FIG. 17. Singlct-gaps A? (sec text) of the iV = 32 M = 4 
sample vs J' for values of J2 « J^{.J')- The dashed lines are 
guide to the eyes. 

three states |2 >, |6 >, |8 > of 5i splits. The gaps A^, 
decrease with increasing J', while Ag increases. The set 
<S2 separate in two groups. The first group consists of the 
states |1 >, |3 >, |4 > which remain very close in energy 
up to J' = 0.6. Their gaps decrease with increasing J'. 
The second group includes the states |5 > and |7 > which 
remain at approximativelly the same energy up to J' ~ 
0.6. Beyond this value, |5 > decreases rather abruptly 
to an energy close to the first group, a drop which may 
be attributed to the proximity of the range of values of 
0.9 < J' < 1 where the crossover to the large J' behavior 
occurs (for J' = 0.9 additional singlet states which evolve 
from the lowest singlet states of L independent chains 
will begin to appear in the low part the spectrum if A'' = 
32). A similar splitting also occurs in most of the sets 
of states which are above ^3 (the states of ^3 remain 
nevertheless quasi-degenerate). The energy of some of 
these states and the energy of the states of iSs decrease 
with increasing J' and become comparable to the energy 
of the the states |5 >, |7 >, |8 > for J' - 0.6. Yet, 
up to J' = 0.4 the evolution vs 1/L of the gaps A? to 
the states of S\ and 52 shows that these splitting are 
much reduced for L ^ 00. All the gaps to the states of 
Si, S2 and ^3 extrapolate to « 0. The singlet spectrum 
thus remains close to the one of independent chains. The 
finite-size results give nevertheless a first indication that 
the splitting of the states in Si, S2 will subsist for L — » 
00. A differentiation in the extrapolated values of the 
gaps becomes more visible at larger J' values. 

Fig. 18 shows the evolution of the singlct-gaps A" vs 
1/L for the M = 4 samples at J' = 0.8 and J2 = 0.45. 
The evolutions of the A^ still remain quasi-linear with 
1/L if L < 8 (iV < 32) on ( J') up to J' = 0.8. A cur- 
vature is hardly visible. But, since a linear extrapolation 



FIG. 18. Singlet-gaps A° (see text) at J' = 0.8 and 
J2 = 0.45 vs 1/L for for iV = 16,24,32 samples of M = 4 
chains The dashed lines are spline fit to the data for iV < 32 
followed by linear extrapolation for A'^ > 32. 

for TV > 32 would leads to negative values for Af, Aj, 
A3, Ag, an upward curvature may be infered for these 
gaps. This affects the precision of the extrapolations for 
L ^ 00 but one may conjecture, in view of Fig. 18, that 
not only A2 still likely vanishes as at J2 = 0.4, but also 
Aj, whereas the other A? are very small and some may 
vanishes. This could be the case of A3, Ag which extrap- 
olate to negative values and probably A4 which follow 
closely A3 for all sizes if J' < 0.8. The vanishing of Ag 
is nonetheless uncertain as the degeneracy of state |6 > 
with the state |2 > seems to be lifted on the line J™ (J'). 
But Aj, Ag and probably A5 which extrapolate above 
A4 remain finite . Thus the degeneracies which occcur 
for decoupled chains are lifted. Nonetheless these gaps 
remain probably very small which suggest that they open 
exponentially like the spin-gap. The gaps to the states 
at the bottom of the singlet continuum are also small. 
Those to the states of .Ss extrapolate to quasi- vanishing 
values. The spectrum appears to be gapless or quasi gap- 
less. 

The exact degree of ground-state degeneracy is diffi- 
cult to acertain, especially for M = 6. The most likely 
vanishing gaps are A?, A§, Ag, Aj^, Ag. for M = 4. This 
would lead to a 2^^^ = 8 ground-state degeneracy. Yet 
the degeneracy would be limited to 6 if Ag remains finite. 
Besides the gaps to the lowest states in the singlet con- 
tinuum appear to be quasi-vanishing. For M = 6, the 
extrapolations of the singlet gaps, although even more 
unaccurate than for M = 4, also indicate a large ground- 
state degeneracy and a gapless continuum of singlet exci- 
tations adjacent to the ground-state. This also suggests 
a spectrum that correspond closer to the one of the RVB 
state of NT than to the one of independent chains. 
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N=32 J' = 0.8 J2=0.45 (a) 




FIG. 19. Dimer-dimer correlation function (see Fig 16) for 
N = 32, J' = 0.80 Ja = 0.45. 

This degeneracy of additional singlet states is accom- 
pagned with a modification of the pattern of the dimer- 
dimer correlations D{i,j;k,l) as shown for A'' = 32 in 
Fig. 19. Dimer-dimer correlations _D(l,5;fc,Z) between 
a reference horizontal pair (1, 5) and a pair (fc, I) on 
the same chain remain strong but have now become ex- 
tremely small otherwise, much weaker than for J2 = 0.4, 
as shown in Fig. 19(a). As shown in Fig. 19(b), this is also 
the case for dimer-dimer correlations D{1, 4; fe, I) between 
a vertical dimer on a pair of chains and vertical dimer 
on a different pair of chains, whereas D{1, 4; k, I) remain 
of significant magnitude for pairs of dimers joigning the 
same chains. Fig. 19(a) (b) give a picture of strongly cor- 
related pairs of adjacent chains, decorrelated from one 
pair to the next with dimer long range correlations along 
individual chains and along the rungs of pairs of adjacent 
chains. This pairing is reminiscent of the RVB state of 
NT. 

The ED results thus indicate a behavior close to the 
one of independent chains on J™ (J') up to very large 
interchain coupling but which differs. Instead, they sug- 
gest that the RVB behavior predicted by NT at weak 
interchain coupling may be realized and also extends at 
large interchain coupling, perhaps up to J' = 0.8, if 
J2 « J^{J'). 

The extension of the RVB state on J™ [J') would how- 
ever have an upper limit. For J' > 0.9 the evolution of 
the singlet gaps vs 1 /L show that they are clearly finite 
if L ^ 00 for Af = 4. This would imply a transition on 
JT{J')- The transition might be discontinouus as the 
gaps seems to open abruptly. An alternative possibility 
could be that the RVB state is only approximativelly real- 
ized on J™(J'). The degeneracy of the RVB state would 
be progressivelly lifted with increasing J' but slower than 
the degeneracy corresponding to the independent chain 



behavior. 

The ED results on the M = L samples of iV = 16, 36 
spins at J' = J indicate that the behavior of the Ji — J2 
model for J™ differs from the one at J' < J. But several 
features reminiscent of this state subsists. The lowest 
triplet excitations have only a small dispersion on the 
boundary of the Brillouin zone. The singlet spectrum 
at finite size is rather similar to the one for J' < J: 
one finds low lying states with wave vector on the side 
of the Brillouin zone in addition to |1 >, |2 >, |3 > 
(degenerate with |2 >) which would allow to form a four- 
fold degenerate columnar VBC. The Ji — J2 model for 
J™, if it displays VBC order, is nonetheless close to a 
state with a spin-liquid behavior at this point of the phase 
diagram. 

C. J2 > Ja™ 

The behavior predicted by NT is also most likely lim- 
ited to the line J^^iJ')- Increasing J2 at fixed J' one 
finds a narrow region with a different behavior before 
reaching the region where the spectra display the fea- 
tures of the (tt, 0) Neel phase. The spin-spin correlations 
s{l,m) (Eq. 5.1) at J' < 0.8 have now the same signs 
as in the (tt, 0) Neel phase: s{l,m) ^ — l', alternating in 
sign in the horizontal direction and being ferromagnetic 
along vertical lines. Fig. 20 shows the evolution of the 
singlet gaps at J' = 0.8, J2 = 0.47 for M = 4 which in- 
dicate that singlet gaps have opened except A" and 
which may still vanish. This excludes an eight-fold de- 
generacy of the ground-state. But the exact degeneracy 
of the ground-state is difficult to acertain. The region is 
very narrow (all singlet gaps are clearly finite at J2 = 0.5 
when J' = 0.8), much narrow than the region between 
the the (tt, tt) Neel phase and the line J2 = J2'{J'). Due 
to the proximity of the line J2 = J2^{J') where A", A" 
would vanish it is difficult to conclude whether both A° 
and A2 really vanishes. 

The patterns of dimer-dimer correlations are shown in 
Fig. 21 for iV = 32 and M = 4 at J' = 0.8, J2 = 0.47. 
The correlations between horizontal dimers in Fig. 21(a) 
display a similar alternation of sign as in Fig. 16(a) and 
could be compatible with a columnar (tt, 0) VBC LRO, 
although they have a small magnitude if the bonds are 
on different chains. 

As the shortest pair of spins with the next largest AF 
spin-spin correlations are on the diagonals of a square, 
whereas the vertical nearest neighbor spins are ferromag- 
netically correlated, we have displayed the correlations 
between diagonal bonds (Fig. 21(b)) and between diago- 
nal and horizontal bonds (Fig. 21(c)). 

As shown in Fig. 21(b), the correlations between diag- 
onal bonds D(1,6;A;,/) do not show a clear modulation 
in their magnitudes neither in the horizontal, the vertical 
or the diagonals directions. Z?(l,6;fc,/) is small except 
for bonds linking two adjacent chains. There subsists a 
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FIG. 20. Same as Fig. 18 at J' = 0.8 and J2 = 0.47. 

tendency of pairs of adjacent chains to associate as seen 
in Fig. 19(b). The correlations are then rather similar to 
those in the 'Haldane phase' of the two-leg ladder. There 
is no visible tendency of the diagonal bonds to order. 

As shown in Fig. 21(c), the correlations between a hor- 
izontal bond and a diagonal bond in the same column are 
negative and there is a sign alternation of the correlations 
between an horizontal bond and diagonal bonds in the 
horizontal direction similar to alternation in Fig. 21(a), 
compatible with a breaking of translational symmetry in 
the horizontal direction. 

A vanishing of A]*, A" would lead to a three- fold 
ground-state degeneracy, but the pattern of dimer-dimer 
correlations do not indicate a VBC LRO with a three-fold 
degeneracy of states |0 >, |1 >, |2 >. 

On the other hand, a columnar VBC LRO with pla- 
quctte modulation, which would lead to an additional 
translational symmetry breaking in the vertical direction, 
appears also unlikely for J' < 0.8 in view of the pattern 
of dimcr-dimcr correlations and the fact that the gap A3 
to the singlet state |3 > with wave-vector (0,7r) is finite. 

If there is VBC LRO, it is most probably a (tt, 0) like 
VBC LRO. It will be associated with a vanishing of A2 
whereas A^ will be finite. The small but noticablc up- 
ward curvature in the evolution of A^ in Fig. 20 supports 
indeed that Aj might be finite for L — > 00 in this region 
for J' = 0.8. As for J2 < J^, the (77,0) VBC LRO dis- 
appears for for J' < 1 if M is finite but could extend for 
J' -» 1 if M = L ^ cx). 

The occurence of a (tt, 0) VBC at J2 < J" (J") would 
mean that there is the same kind of VBC LRO on both 
sides of J^i^J') for M > 4. The state for J2 > J^i^J') 
could however differs from the state for J2 < J^{J'^ 
by some non-local order parameter similar to those of 
the 'singlet state' and the 'Haldane state' of the two-leg 
ladder as conjectured at the end of Sec. IV for the phases 



N=32 J' = 0.8 J2=0.47 (a) 




FIG. 21. Dimer-dimer correlation functions (see Fig 16) for 
AT = 32, J' = 0.80 J2 = 0.47. 

at finite M > 2 on both sides of the intermediate region. 

VI. SUMMARY 

Exact diagonalization calculations provide support 
that the RVB behavior predicted by Nersesyan and Tsve- 
lik in the J —J' — Ji model at J2 = 0.5 J' <C J may be re- 
alized and extends along a curve coincident with the line 
of maximum frustration J™ (J') where it subsists up to 
large interchain coupling for an even number M of chains 
of lenght L ^ 00. The line J2'{J') is the analogue of the 
transition line J2 = 0.5 J' in the classical model on which 
the chains decouples. The independent chain behavior 
appears destabilized by quantum fluctuations. 

The line J 2^ {J') is located in an intermediate region 
of the phase diagram where the present results suggest 
that the ground-state displays a columnar (tt, 0) valence 
bond long range order. This region of columnar order 
will have a finite width as soon as M > 4 and extends on 
both sides of J™ (J') up to a value of J' close but smaller 
than J for finite M, that increase with M and approach 
J' = J for M = L ^ 00. It separates two phases that 
are fully gapped Haldane like phases for finite M and 
evolves respectivelly to the (tt, tt) and (tt, 0) Neel phases 
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for M ^ oo. The Haldanc like phases have probably a 
topological order. Its nature for M > 2 deserves further 
investigation. 

The occurence of a columnar order for J' smaller than 
J, next to the Neel phase, agrees with the scenario de- 
duced from large- approaches that valence bond crystal 
order is expected after the destabilization of a collinear 
AF Nccl phase. It is possible that the columnar order 
may extend for J' ^ J if M = L ^ oc and that this 
scenario may be realized also in the Ji — J2 model which 
corresponds to J' = J. The behavior of the J — J' — J2 
model in the intermediate region for J' ^ J requires 
nonetheless further study to confirm this hypothesis. 

The exact diagonalization method is limited to systems 
of small sizes. The extrapolations to the infinite lenght 
limit are affected by uncertainities. Indeed, some one di- 
mensional systems requires data at very large lenghts to 
allow a reliable extrapolation to the infinite lenght limit. 
For certain systems there may be a crossover in the scal- 
ing behavior of the gaps at a very large lenght as shown 
for instance in Ref. [38] where some gaps do not decrease 
monotonicaly with the lenght beyond a certain lenght. 
For the present model we have assumed a monotonicaly 
decreasing smooth evolution of the gaps with increasing 
lenght. The assumption appears reasonable in view of the 
data available at large lenghts in the case of two chains 
(the 2-leg ladder). Yet, a small change in the evolution 
of the gaps at very large lengths is not to be excluded 
for those gaps that appear to extrapolate to vanishing 
or quasi vanishing values. The scaling behavior of the 
gaps to the low lying excited states as a function of the 
number of chains deserves further studies. In particular, 
calculations using an approach that can deal with very 
large systems, like possibly the DMRG method, would 
be worth to firmly establish wether the gaps that have 
been conjectured to extrapolate to zero really vanish or 
are just small but finite in the thermodynamic limit and 
thus check the exact degeneracy of the ground state on 
the line of maximum frustration. 

This study was limited to the case of AF couplings. 
Yet, the J — J' — J2 model with AF intrachain coupling 
(J > 0), but interchain couplings (J', J2) that arc either 
or both ferromagnetic (FM) is also of interest [39]. In- 
deed, FM couplings may be relevant to modelize certain 
2-lcg ladder materials [40] and might be introduced to 
describe some M-leg compounds. 

Classically, the ground-state of the J — J' — J2 model 
remains AF for J' > -1 and J2 > -0.5 (if J = 1). It 
is the (7r,7r) Neel state for J2 < 0.5 J', the (7r,0) Neel 
state for J2 > 0.5 J', and the J chains stay decoupled for 
J2 = 0.5J'. 

The spin-1/2, 2-leg ladder (with FM J' and/or J2 ) 
has been studied by bosonization [41] and DRMG [42] 
calculations. These studies showed that the singlet and 
Haldane phases extends for J', J2 < 0, separated by a 
transition line which is the continuation of the bound- 
ary line for J', J2 > 0. This transition line is also the 
line of the points of max;imum frustration J™ {J'). The 



transition line, however, does not deviate from the clas- 
sical boimdary line J2 = 0.5J' and the transition ap- 
pears to be 2nd order for J', J2 < 0, instead of 1th order 
for J', J2 > 0. In the bosonization approach, the inter- 
chain interaction, derived at 1th order, proportional to 
g = J' + 2J2 which now turn to be < 0, is then irrelevant 
and does not open a gap. The system is predicted to be 
critical and to consist of independent chains [43] . 

A similar change of sign of the interchain interac- 
tions appears in the bosonization approach of NT if 
J2 = 0.5J' < for the 4-leg ladder. The same critical 
behavior could occur on the transition line in the 4-leg 
ladder and perhaps for M-leg ladders {M > 2 even) . Pre- 
liminary calculations indicate that the curve J^^^J') also 
remains on the classical boundary line for M = 4, where 
the model might be critical, perhaps down to J' = —0.8. 
The model may exhibit "sliding Luttinger spin-liquid" 
behavior as proposed for the crossed chain model [44], 
at least for not too strong interchain interactions. On 
the line J2 = 0.5J', the lowest energies E{S) in spin 
sector S evolves linearly with S. For M ^ 00, there 
may be also an intermediate region with a spin-liquid 
state around this line in between the Neel phases, where 
E{S) — Eq ^ S{S+\). This issues deserve further study. 
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